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AVERTISSEMENT, ' ^ 

Vi 


Je me suis attaché k ne donner à ce Complément 
que le peu d’étendue strictement nécessaire à l’objet 
que j’avais en vue; seulement, quoique le Pro- 
gramme se borne à exiger la résolution des équations 
du troisième degré, j’ai préféré appliquer la méthode 
des substitutions par les différences à la séparation 
des racines d’une équation du cinquième degré, sans 
que le nombre des racines réelles pût être connu à 
priori, afin de mieux faire apprécier cette méthode, 
et de montrer qu’elle ne laisse subsister, en général, 
aucune incertitude. J’ai effectué le calcul des valeurs 
approchées des racines par des procédés simples et 
rapides. Il est d’ailleurs facile de juger du degré de 
l’approximation qui a été obtenue, sans recourir à 
des règles générales dont les prescriptions absolues, 
loin d’offrir aucun avantage réel, ne seraient souvent 
qu’une gêne pour le calculateur. Je pense que les 
exemples que j’ai traités paraîtront suffisants pour 
mettre les élèves à même d’acquérir la pratique qu’on 
leur demande, relativement aux équations algébri- 
ques et transcendantes. 

. Avril i85i. 
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COMPLÉMENT 

D’ALGÈBRE. 


Formules générales relatives aux différences . 

1 . Considérons une suite de quantités 

Uof W|f U3 y . • • • 

Si l’on retranche chacune de ces quantités de celle qui la 
suit immédiatement, on formera une autre suite dont les 
termes seront ce que l’on appelle les différences premières 
de la suite proposée. En désignant ces différences par la 
lettre A placée devant la quantité que l’on soustrait, on a 

au„ = m, 

= «, 

A U, = «3 


Si l’on forme de la même manière les différences des 
quantités Au,, Au,, Au,, etc., on obtiendra les différences 
secondes des quantités proposées, qu’on désigne par A*. 
<3n a ainsi 

A’tt B = Au, — Au„,. 

A’u, = a«, — Au,, 

A’ u, = A u 3 — A u„ 


— «0, 

— «„ 
— «„ 


En continuant de même, on formera successivement les 
différences troisièmes, les différences quatrièmes , etc., qui 
seront désignées par A 3 , A*, etc. 
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COMPLÉMENT D'ALGÈBRE. 

2. On peut obtenir facilement l’expression d’une diffé- 
rence d’un ordre quelconque , en fonction de toutes les 
quantités de la suite primitive dont elle dépend. On a d’a- 
bord, par les valeurs ci-dessus de A « 0 et de Au,, 

. A 5 M 0 = II, — U I — («, — Mo) = U, — 211,+ M,; 

pareillement 

A 1 /!, = Wj — ?.«, + u, : 

donc 

A 1 !!, = («j — 2U, -+-«,) — (il, 2M, -+- II,) 

= M, — 3 U : — |— 3 II ! — Mo- 

Les coefficients numériques de ees expressions de A* u n 
et A n u 0 sont les mômes que ceux du carré et du cube d’un 
binôme dont le second terme est négatif, ce qui donne lieu 
de penser que l’on aura généralement 


. , n in — i) 

( I ) A = il,, — -t- — M„ 


.2 


n{n — Q(« — 2) . 
1.2.3 


On peut, en eflét, démontrer cette formule, en formant 
l’expression de A n u 0 , au moyen de celles de A"~ 1 u 0 et de 
A"' 1 u, . Soit , 

• * -/ / 

A"- 1 u, = ii« — i + B u n — j — C 4 -+- , . . ; 


on aura pareillement, et avec les mêmes coefficients, 


A"-’ m, = m„ — A ti n . ., -t- B tifx-Li — Cm„_, -t- . . . , 


et de là on conclura 

A" Mo — A"-* 1 U, 'A" -1 Mo 


II 

à 

1 

> 

ii« — i B 

m„_ : — G 

1 

-t-A 

1 — B 


On voit que les coefficients numériques , dans le passage 
d’une différence à la suivante, se forment suivant la même 
loi que dans le passage d'une puissance du binôme x — u à 
la puissance .suivante-, et puisque les coefficients numé- 
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COMPLÉMENT d’aLCÈBRK. 

tiques do l’expression de A* u 0 sont les mêmes que ceux du 
carré de x — a , les coefficients de l'expression d’une diffé- 
rence quelconque seront aussi les mêmes que ceux de la 
puissance correspondante. 

La formule (i) peut être écrite de cette manière symbo- 
lique , 

A"it = (rt — i)", 

sous la condition qu’en développant la puissance indiquée, 
on devra changer les exposants de u en indices. 

3 . Lorsque l’on connaît les différences de l’ordre tt d'une 
suite de quantités, et le premier terinc de chacune des 
suites précédentes de différences, et de la suite des quantités 
rHes-mèmes , on peut former toutes ces suites. On obtien- 
dra d'abord celle des différences’ de l’ordre n — i en ajou- 
tant successivement à la première de ces différences, qui est 
connue, celles de l'ordre n ; on passera de même de la suite 
des différences de l’ordre n — i à celle des différences de 
l’ordre n — 2 ; ainsi de suite. 

A. On peut aussi obtenir l’expression de l’une quel- 
conque des quantités de la suite primitive, au moyen du 
premier terme u 0 de cette suite et de ses différences succes- 
sives An,, A* n*, A n i/ 0 , etc. 

On a » . 

«, = il, -4- a Hj «t Au, — Ait 0 + A’ U , ; 

H, =r «i -+- A U, = U, 4- 2 A U, + A 5 U , , 

A 11,=: A II, -4- A 2 «, = A u, -f- 2 A’ il, -+- A J «„ , 

«j = «j - 4 - Au, = u, -t- 3 Ait,, 3 A 1 11, -+- A 

Les coefficients numériques des expressions de «, et de n, 
sont encore ceux du carré et du cube d’un binôme, cl il 
est facile de démontrer que celte loi est générale. Si l'on 
suppose 

Hfi—, - u, -t - A Au, -4- 11 A’ /t 0-4- C A 1 «0 - 4 - . 1 . , 
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4 COMPl.KMSNT D AK.KBRP. 

on aura pareillement 

Aa„_i = Au, 4 - Ai 1 #, 4 - BA'k 0 4 - CA'«, 4-. . . , 


et il en résultera 

u n = i 4- Au„_ 


= U, -t- A 

Au, 4 - B 

A 5 «0 C 

•4” 1 

4- A. 

+ B 


A 3 w» 4~ . , 


en sorte que les coefficients numériques, dans le passage 
d’une quantité à la suivante, se forment de la même ma- 
nière que dans le passage d’une puissance de x -4- a à la 
puissance suivante. On aura donc généralement 


(a) b 1 + ijAb,,4 




A' u, - 


n(n — i ) ( « — a) 
i . a. 3 


A «„ 


Cette formule peut être écrite symboliquement de cette 
manière , 

«„ = (i 4- A) n u; 

il faudra seulement observer qu’ après le développement 
de la puissance, les exposants n’auront pas leur signification 
ordinaire , et ne seront que des caractéristiques de l’ordre 
des différences. * . ‘ 


Différences des fonctions entières. 

5. Supposons que les quantités u 0 , «j, Mj, « 3 , etc., soient 
les valeurs d’une fonction entière de a:, et qu’elles corres- 
pondent à des valeurs de la variable X , équidistantes les 
unes des autres. Si la fonction est du degré n , elle sera 

Af 4 Bî *" 1 4- C.r "- 2 4- D 4- 

Désignons par h la différence de deux valeurs consécu- 
tives de Xj on aura, en représentant par u la fonction , 

Au = A [.(.r - 4 - A ) " — •*"] 4 B [ (.t 4 h) n ~' — x”~' 
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COMPLEMENT 1» ALGÈBRE. 5 

Cette expression Je Au est du degré n — i par rapport à x , 
et de cette forme , 

nAhx"-' B, hx n ~' 1 4- C, Ax' 1-3 4- . . . . 

On obtiendra la différence suivante du même ordre, A 
en remplaçant x par x ■+■ h ; et , puisque A s u = Au, — Au, 
on trouvera 

A’ u — n (n — i) A A’** - ’ -t- B, h 1 x" -3 4- C 3 A 1 x" _ ‘ 4- . . . . 

Les degrés de ces différences successives Au, A * u, etc., 
décroissant d’une Unité, de chaque différence à la suivante, 
la différen'ce de l’ordre n ou A" u ne dépendra pas de x , et 
elle sera 

X‘ u — i . a . 3 . . ..n A h n > 

Cette différence étant constante, toutes celles des ordres 
suivants seront milles. 

G. D’après cette propriété des fonctions entières , il est 
facile d’obtenir la Suite des valeurs d’une telle fonction , 
pour des valeurs équidistantes de la variable , quand on a 
calculé directement un nombre de termes de cette suite égal 
au degré de la fonction. 

Soit-, par exemple, le polynôme 

x s — 3x 3 4-a; î — 8x — io; 

et supposons qu’on veuille obtenir les valeurs de ce poly- 
nôme en prenant pour x les nombres entiers positifs et 
négatifs. 

Il faudra calculer directement les valeurs correspondantes 
à cinq nombres entiers consécutifs , et l’on pourra choisir, 
pour plus de simplicité, les nombres — 2, — i ,0,4- i, 4- 2. 
On formera ensuite, avec les résultats de ces substitutions , 
le tableau suivant : 
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• La première colonne' contient les valeurs de x , et la 
deuxième, les valeurs correspondantes de la fonction.- On a 
d’abord formé les nombres qui sont, au-dessus des filets , et 
l’on a obtenu ceux de la troisième colopne en retranchant 
chacun des nombres de la deuxième colonne de celui qui 
est au-dessous de lui. La quatrième colonne est formée avec 
la troisième, comme celle-ci a été formée avec la deuxième. 
La difierence cinquième est connue à priori par ce qu ? on 
a vu dans le numéro précédent. 

Pour obtenir les nombres qui sont au-dessous des filets , 
on continue d’abord la colonne des différences quatrièmes , 
dont chacune est égale à la précédente augmentée de la dif- 
férence constante 120 . On forme chacune des dilïërences 
troisièmes en ajoutant à celle qui précède la différence qua- 
trième placée dans la môme ligne. On agit de la même ma- 
nière pour obtenir les différences secondes, les différences 
premières, et enfin les nombres contenus dans la deuxième 
colonne à gauche, qui sont les valeurs de la fonction. 

11 faut opérer autrement pour les valeurs correspondantes 
aux nombres négatifs — 3, -r-4? etc., qui sont données 
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par le tableau ci-dessous : 


■r = — 3 

- «fy l 


141 

- 14 * 

-i- l32 

— 130 

— 4 

- 79'» 


655 

— 5>4 

-+- 372 

— 24° 

— 5 

— 2695 

-h 

njoi 

— 1246 

-t- 732 

— 3(io 


Ou retranche d’abord la différence constante 120 de la 
différence quatrième o correspondante à x — — 2 , ce qui 
donne la différence quatrième — 120, correspondante à 
x = — 3 . On retranche cette différence quatrième — 120 
de la différence troisième 4- 1 2, correspondante à x = — 2; 
ce qui donne la différence troisième + i 32 , correspondante 
à x = — 3 . On obtient d’une manière semblable les diffé- 
rences — i 4 2 î ■+" t 4 i et le résultat — i 3 g. Les nombres 
de la ligne correspondante à x = — 4 son! formés avec 
ceux de la ligne correspondante à x = — 3 , comme ceux- 
ci ont été formés avec les nombres contenus dans la ligne 
correspondante k x — — 2; et l’on continue de même autant- 
qu’il est nécessaire. 

7 . Lorsque deux fonctions entières, du degré n , ont les 
mêmes valeurs pour n H- 1 valeurs différentes de la va- 
riable, elles sont identiques; car, si elles ne l’étaient pas, 
leur différence , qui serait du degré n au plus , serait an- 
nulée pour n 4- 1 valeurs de x : ce qui est contraire à la 
proposition qu’une équation du degré n n’a pas plus de n 
racines. . _ . 

En représentant la fonction par 

À*» +Rj"-' + C x "- 1 

et en exprimant qu elle devient successivement h 0 , u , , 
t< 3 , etc., lorsqu'on remplace x par les valeurs désignées x 0 . 
Xi, x t , etc., on aurait n 4 - 1 équations , au moyen desquelles 
on trouverait les valeurs des n 4- 1 coefficients A , B, C, etc . 

Mais , dans le cas où les quantités x 0 , x,, x, . . . , x„ 
sont équidistantes les unes des autres , l’expression de la 
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fonction peut se déduire de la formule (a), [n ü 4j. Le se- 
cond membre de cette formule devenant successivement h 0 > 
a., «j, etc., lorsqu’on fait n = o, = i , = 2 , etc. , si l’on 
représente par h la différence constante des deux termes con- 
sécutifs de la suite x 0 , Xi , .r., etc., et si l’on écrifdans la for- 
mule ( 2 ), au lieu de n, - - 1 le second membre deviendra 

successivement i/ 0 , « s , . . . , u„ , lorsqu’on rempla- 

cera x par .r 0 , x 0 h, x 0 4- a/», . . . , x„ -+- tili. D’ailleurs 
il sera une fonction entière de x , du degré «; il sera 
donc la fonction cherchée. Ainsi , en désignant cette fonc- 
tion par u x , on aura 


(3) 


u z = u, -f- 



A u 0 -f- 




A* u p 
I . 2 



A J //„ 

ô “H 

1.2.3 


8. Dans le cas particulier où l’on fait x = x 0 4- J A) on 
obtient pour la valeur de la fonction, 

(4) «„ + jâu,- j a 1 «, -f- -Jj A 3 a, — a 1 «, 4- 


9. Si l’on fait, pour abréger, - — -j~ — z , la formule (3) 
devient 

z(z — i) z(z — i)(z — 2 ) 

(5) U x — u a -*-zXu„ 4-— V H — — ■ A 3 M„.4- .... 

1.2 1.2.3 


Supposons qu’on ne donne à z que des valeurs comprises 
entre o et 1 , ce qui revient à donner à x, dans la formule (3), 
des valeurs comprises entre x 0 et x„ 4- h. Dans ce cas, les 
coefficients des différences Ai/ 0 , A' « 0 , A 3 h 0 , dans la for- 
mule (5 ) , sont alternativement positifs et négatifs , et leurs 
valeurs absolues vont eu décroissant. Le coefficient de 

A % — Z ^ ^ > a une valeur au plus égale à A ; car, la 

somme des deux facteurs z et 1 — z étant constante et égale 
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* 9 

à i , le maximum de leur produit est j. Ou détermine aisé- 
ment le maximum du coefficient de A 3 , et l’on trouve qu’il 
est moindre que Il en résulte que la valeur numérique 
du coefficient de A 1 est moindre que ^ X f ou ; celle du 
coefficient de A 5 est moindre que ^ X { ou etc. 

10 . Je vais donner un exemple de l’application de ces 
principes à la séparation des racines d’une équation numé- 
rique. Soit l’équation 

x i — 3 .r 5 x 1 — 8 x — io = o. 

Les résultats de la substitution des nombres entiers, qui 
ont été formés précédemment (n" G ) , fout connaître qu'il y 
a au moins une racine entre les nombres 2 et 3 ; qu’il s’en 
trouve une entre o et — i , et une entre — 2 et — 3 . On peut 
aussi s’assurer qu’il n’en existe aucune en dehors de ces 
trois intervalles. Comme on voit, par la règle des signes de 
Descartes, que l’équation n’a pas plus de deux racines néga- 
tives, il suffit de considérer les résultats de la substitution 
des nombres positifs. 

Pour l’intervalle des nombres o et i, en donnant dans la 
formule ( 5 ) à u„ la valeur correspondante à x — o, et qui 
est — îo, les termes positifs ne proviendront que des diffé- 
rences troisième et cinquième, égales à a 3 a et 120; mais 
011 juge aisément qu’en raison du terme dépendant de la 
dillérence quatrième, 240, qui sera négatif, la fonction u x 
11e pourra pas changer de signe (*). 


(*) On peut cci-ire ainsi les trois derniers termes de la formule (5) : 


j (1 — (a — j) 




t (l— t)('J — z \ (3 — z\ — z) 


z clam moindre que 1 , d’api ès 

1 — 3 

y— A* obi moindre qUc 


les valeurs A* = i3a, A*'= ‘i\ o , la quaniilc 
ri, ei, par conséquent, la première partie 
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Par rapport à l’intervalle des nombres i et i , ou recon- 
nail encore que les trois derniers termes de la formule (5) ne * 
peuvent donner une quantité positive égale à la valeur nu- 
mérique de n 0 q u * csl — 19 ; à la vérité, la différence pre- 
mière donne un terme positif ; mais , pour que ce terme 
positif puisse l’emporter sur celui qui dépend de la diffé- 
rence deuxième , et qui est négatif, il faut que la valeur de x 
diffère très-peu de l’unité , et cette condition rend très- 
petits les coefficients des autres termes positifs. 

Pour l’intervalle des nombres 3 et 4i la formule (5) 
donne 

iij = 1 374-3 [ 66 g — (t — 3)625] 4 -z(i — z)(2 — z)[202 — ( 3 — 3)26] 

4-3 (r — 3)'(2 — 3 ) (3 — 3 ) (4— 3 ); 

z étant moindre que 1 , tous les termes de cette expression 
de u x . sont positifs. 

O 11 reconnaît de la même manière qu’il n y a pas de ra- 
cine comprise entre 4 el 5. 

Au delà de 5, en prenant pour la valeur t/ 0 celle qui 
répond à x = 1 , on a 3 = 4 ou 3 )> ] ; tous les coefficients 
des différences sont positifs", et, d’après les valeurs de ces 
différences qui correspondent à x = 1 5 on voit que la somme 
des termes positifs l’emportera constamment sur le premier 
terme, qui sera seul négatif . 

11. Il est mile de connaître les valeurs des coeffi- 
cients des différences, dans la formule (5), pour le cas 
où l’on donne à z les valeurs 0 , 1 ; 0 , 2 ; o,3 ; . . . , 0 , 9 . 
Ces coefficients étant des fonctions entières de z , on peut 
former leurs valeurs au moyen des différences successives; 


de la quantité ci-dessus est moindre que la seconde partie est plus pe- 
tao 

lt le nue — — • 

1 M'I 
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elles sont telles qu’on les voit dans le tableau ci-dessous : 


A 

à 3 

A 3 

A 1 

A 1 

-h 



“h 


-h 

0,1 

0,045 

0,0285 

0,020.66a . 5 

o-,oi6. 116.70 

0,2 

0,080 

0,0480 

o,o 33 .Goo.o 

0,025.536.00 

o ,3 

0,1 o 5 

o,o 5 i )5 

0,040. iG*j. 5 

0,02<J.720.!î5 

0,4 

0,120 

0,0640 

0,041 .Goo.o 

0,029.902.0 0 

o,5 

0 , 1^5 

0,0626 • 

0,039.062.5 

0,027.343.70 

o,G 

0,120 

0 , o. r >Go 

o,o 33 . Goo.o 

0,022.848.00 

0.7 

0,1 o 5 

o,o455 

0,026. 162.5 

0,017.267.2s 

0,8 

0,080 

o,o 3 ?o 

0,017. Goo.o 

0,01 f .26.4.00 

0,9 

0,045 

o,oi 65 

0,008.662.6 

0,005.370.7.5 


Au moyen des différences qui ont servi à calculer ce 
tableau, on obtient immédiatement les différences de la 
fonction u x pour des valeurs de s, de dixième en dixième, 
et qui correspondent à des valeurs de X croissantes par difié- 

rences égales à ~ En désignant ces nouvelles différences 

par la lettre d, on trouve : 


â’u 0 = 
S‘u,= 

$>u,= 


A — o,o 45 A ! -f- 0,0285 A 3 — 0,020662.5 A' 4 - 0,0161 i 6 .^ 5 a 5 -h.. 
o, 01. A 3 — o,oot).A 3 -t- o, 0077254' — 0,006697 . 5 a s +... 

0,00 1 A 3 — o,ooi 35 a' -h 0,001462.5 A 3 4 -... 

, 0,0001 A* — o,oooi8a s 4-.... 

0,00001 A 4 4-... 


Ces formules et le tableau qui précède suffiront jusqu’au 
cinquième degré inclusivement. 


12 . En appliquant les dernières formules à l’exemple du 
n° 10 , pour substituer Iesdifférentsnombresdedixièmesentre 
2 et 3 ,on trouvera d s « 0 = 4 - 4 , 668 . 01, d ! // 0 = -(- 1,496.3, 
d 3 /t„= 4-0,259.5, d*«« = 4- 0,026’. 4, cî 5 = 4-0,001.2. 
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On formera ensuite , avec ces différences , le tableau 
suivant : 


x =: 2 

— 14,000.00 

-+- -4.668.oi 

-t- 1,496.3 

-+- 0 , 2 J 9.5 

-H 0,026.4 

-t- 0,001.2 

2,1 

— 9. 33, -99 

-t- 6 , 1 64 . 3 1 

H- 1 , 755.8 

-t- 0,285.9 

-1- 0,027.6 


V 

- 3,167.68 

-t- 7,920.11 

■+- 2,0.41 .7 

-+- o, 3 i 3.5 

-+- 0,028.8 


2,3 

-t- 4 . 752.43 

-t- 9,961.81 

H- 2 , 355.2 

0,342.3 

-t- o,o 3 o.o 


a ,4 

-+- i4.7«4-a4 

-H 12, 317. 01 

-t- 2,697.5 

-t- 0,372.3 

-i- o,o3i.2 


•j , 5 

-+- 27,o3i.i5 

-H i 5 ,oi/|. 5 i 

-h 3 , 069 8 

-t- o, 4 o 3.5 

-+- o,o32-4 

• J 

2 ,6 

-+- 42,055.76 

-t- i 8 ,o 84 . 3 i 

H- 3 , 473.3 

H- 0,435.9 

-t- o,o 33.6 


a» 7 

-H 60,130.07 

-1- 21 ,557.61 

■+■ 3 , 909-2 

•+• 0,469.5 



a, 8 

-t- 81,687.68 

-1- > 5 , 466 . 81 

-t- 4,378.7 




a .9 

-+- 107 , 154.49 

-1- 29 , 845.51 





x =. 3 ,o 

“t" 137 , 000.00 







On conclut de ces résultats qu’il n’existe aucune racine 
en dehors de l’intervalle des nombres 2,2 et 2,3; car, 
dans chacun des autres intervalles, les termes de la for- 
mule (5) qui auront un signe contraire à celui du premier 
terme, ne pourront donner leur signe à la valeur de la 
fonction. 

On verra , dans les articles suivants, comment 011 peut 
parvenir, par des moyens plus rapides, à des valeurs plus 
approchées de la racine comprise entre 2,2 et 2,3. 

Si l’on avait eu la certitude que l’équation n’admettait 
pas plus d'une racine positive, on aurait pu arrêter le ta- 
bleau ci-dessus au résultat de la substitution de 2,3. Mais 
alors, on n’aurait eu aucune garantie de l'exactitude des 
opérations, tandis qu’eu les prolongeant jusqu’à la valeur 
correspondante à a; =3, on obtient une vérilication, 
puisque cette valeur est déjà connue. 

Exemple du calcul de la -valeur approchée d’une racine 
d’une équation numérique. 

13. En développant le second membre de la formule (5), 
et en faisant «,. = «>, on obtient une équation de cette 


Digitized by Google 


COMPLÉMENT Ii’aI.GKBRR. 


forme : 

( 6 ) o u ~t- A z — t- B z 3 — (— C z* -f- P z ! E z ^ -4- .... 

On a d’ailleurs 

A=i-{i ! + i4 1 -Ti l +-r i> >" 

B = 7 A 1 • — yi’ + lji 1 — . . 

c = 7 A 3 i-A«r+-^A s ... 

T) — -!- A' !- A‘ .*. . 

E= -J— A 5 . . . 

Lj u #* 4 ••• 

On formera aisément tous les termes de ces dernières for- 
mules , suivant le nombre des différences qu’on aura à consi- 
dérer, et qui dépendra du degré de l’équation proposée. Les 
expressions ci-dessus comprennent tous les termes qu’on aura 
à calculer, dans le cas d’une équation du cinquième degré. 

En négligeant, dans l’équation ( 6 ), les puissances de z 
supérieures à la première, on aura 



Si l’on néglige seulement les puissances de z supérieures à 
la seconde, on aura la valeur de z par une équation du se- 
cond degré. 

1 i. Puisqu’on a fait X ■ ■■ X ~ == s , d’où x = .r 0 -t- liz , si 

l’on représente l’équation proposée par f (x) = o , l’équa- 
tion ( 6 ) sera la même que celle-ci : 

/(*„) -h/' (x a )hz-tf" (.r„) o. 

La correction donnée par la formule ( 7 ) est donc celle que 
fournit la méthode de Newton ; à celte seule différence près , 
que la valeur de z devra être multipliée par h. Mais, au 
lieu d’employer les dérivées du premier membre de l’équa- 
tion proposée pour former l’équation en z, il est plus com- 
mode de se servir des différences qui ont été déjà calculées. 

15. Dans l’exemple qui a servi plus haut , on a , pour la 
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racine positive, qui est comprise entre 2,2 et 2,3 : 

«= — 3,167.68, A=6,9Ç)6.8, 8 = 0,876.8, 

C = o,o 45 . 4 1 D = 0,001.1, E = 0,000.01. 

La valeur de z donnée par la formule ( 7 ) est trop grande, 
puisqu’on néglige des termes qui sont tous positifs. Cette 
valeur trop grande étant moindre que -j , la somme des ter- 
mes négligés' est moindre que 7 B -+- \ C ■-+- -A D -H- ■— E. 
Cette quantité est moindre que 0 , 225 ; ainsi l’erreur de la 
valeur de z est moindre que ; elle est donc moindre 

queo,o 4 , et l'erreur de la correction o,iz est moindre 
que 0,004. - 

J6. Calculons la valeur de z d’après les trois premiers 
termes de l’équation (6) qui donnent une équation du se- 
cond degré. On ne doit prendre que la racine positive; et , 
en changeant u en — u, attendu que la valeur de u est néga- 
tive, on a • • 


En effectuant le calcul au moyen des logarithmes, on 


trouve 

z = 0,429604. 

Voici le tableau des opérations : 


log « 

o, 5 oo. 74'3 


log B 

1 ,942.9005 


B« 

o, 443 . 64 i 8 

2,777.42. 

l»g 7 A 

0,543.8695 


(7 AO’ 

1,087.7390 

12,238. 80 

(i A)* + B«. 

1,176.5607 

l 5 , 0(6.22 

V 

o,588.28o3 

3,875 077 

7A.. 


3,498.4 


1,575.9691 
1,942 qoo 5 

0,376.677 

log - — 

1 , 633 . 0686 
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17 . Si l’on veut pousser plus loin l’approximation, on 
fera dans l’équation (6) z — y -t- z', en désignant par y la 
valeur de z déterminée par la formule (8). On obtiendra 
l'équation 

(9) u' -H AV -H BV’-HÇ' i'’ + DV‘ -t- E' z' 1 = o; 

et l’on aura 

«' = C y 3 H- D 7' -H E y*, 

A'= A-H 2B7 4 - 3C7 1 -H 4 D 7 : '-H5E7‘, 

B f “ B -H 3 C 7 -f- 6D7’ -H 10 Ry*, 

C' ™ C -H 4D7 -H 1 o E7 ? , 

D'= D-H 5 Ev, 

E' = E. ; 

En négligeant toutes les puissances des' supérieures à la 


première, on obtient 



On reconnaît à l’inspection des expressions ci-dessus de n'et 
de x¥, et des valeurs des nombres A, B, C , D , E, y, que la 
quantité z' déterminée par l’équation (10) sera moindre 
que o,oot ; et comme ou ne peut pas être assuré, en opé- 
rant par les logarithmes, d’avoir la valeur de 7, qui est 
celle qui a été calculée dans le numéro précédent , avec plus 
de six chiffres exacts , il serait superflu de calculer la quan- 
tité z! avec plus de trois ou quatre chiffres ; en conséquence, 
on pourra abréger le calcul de u! et de A', en ne prenant 
les logarithmes qu’avec cinq chiffres décimaux seulement : 
on fera donc ce calcul comme il est indiqué ci-après : 
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Calcul de u . 


Calcul de A'. 


logC... 2, 65^.05 A 6,9968 

y '... 2,899.20 2B7. .. 0 , 753.35 

3,556.25 0,003.599.6 Iog3C.. i,i34.i8 

logD... 3,041.39 y 7 . .. I,266.l3 

y '... 2,532.27 2,4oo.3i o,oî5. 14 

5,673.66 0,000.037 5 log4Dy’. 4»54 2 -65 o,ooo.35 

log E y s .. 7,165.34 0,000.000.1 7,775.64 


u' 

logA'.. . 


3,560.77 0,003.637.2 
0,890.74 


4,670.03 


-, = 0,000.467.8 


O11 a négligé, dans le calcul de A', le terme 5Ey s , attendu 
qu’il ne peut influer sur les chiffres dont on tient compte 
dans les autres termes. 

On conclut du dernier résultat et de la valeur de y qui a 
été trouvée précédemment, x = 2,242.91.3.62. 


18 . On peut apprécier comme il suit le degré d’exacti- 
tude de ces approximations successives. Pour la première 
correction , que nous avons désignée par y, on obtient une 
quantité trop grande. Mais les termes de l’équation (6)qu’on 
a négligés étant ceux qui donnent, dans le calcul de la cor- 
rection suivante ,1a quantité u 1 , qui est moindre que 0,0037, 
la valeur de s donnée par la formule (8) deviendrait trop 
petite, si l’on diminuait la valeur de u de 0,0037. Or, en 
reprenant le calcul de z avec la valeur de u ainsi diminuée, 
on trouve que celle de z surpasse 0,42917. L’erreur de la 
valeur de z, qui a été calculée en premier lieu, est donc 
moindre queo,ooo 5 . 

Par rapport à la seconde correction, la valeur de z’ étant 
négative, on peut changer z' en — z' dans l’équation (9). 
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qui devient alors 

(n) aV-4-z'’(B'— CV)- t-z'‘(D'— EV) = o. * 

Les quantités B' — C'a' et D' — E ' z' sont positives; on 
outre , la quantité y étant moindre que j , on a 

B’<B-+4C + iD + {E, D'<D + |E; 

donc 

B'<^o,<)466, D'<^o,ooi4- 


r • 

La correction —, < 


f M 1 


quation (n) que l’on néglige sont positifs; mais la valeur 
de z' devant être plus petite que o,ooo 5 , l’erreur est 


moindre que • 


tité ne peut pas influer sur les six premiers chiffres déci- 
maux de la valeur de z' . On est donc en droit de compter 
sur l’exactitude des sept premiers chiffres décimaux de la 
valeur de x conclue des deux corrections. 


Si l’on prend pour z le nombre 0,429136, en mettant 
en même temps dans l’équation (6) les valeurs de u , A , B, 
G, I), E (n° lo), et en calculant les termes qui dépendent 
de z avec le secours des logarithmes, et jusqu’au sixième 
chiffre décimal, on trouve que la somme de ces termes 
est 3 , 167673-, ainsi l’on a un résultat négatif : mais en pre- 
nant z = o, 429137, le logarithme de z augmente de iouni- 
tés du dernier ordre, et la somme des deux termes en z et z* 
est augmentée par là de o , 000.007.4 > de sorte que l’on a un 
résultat positif. On est assuré par cette vérification qu’il n’a 
pas été commis d’erreur dans les calculs, et que la valeur 
de z est exacte jusqu’au sixième chiffre. 


19 . Il faut remarquer qu’après le calcul de la première 
correction tel qu’il a été exécuté (n° IG) , il aurait été super- 
flu de calculer la valeur de z' au moyen des trois premiers 
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termes de l’équation (y); car, à cause de la petitesse des 
nombres u' et 11', la valeur de z', calculée de cette manière, 

dillererait peu de — — >> et la différence serait inappréciable 

au degré d’approximation de la valeur de z précédemment 
calculée '(*). 

20. On voit que dans le calcul de la première correction y 
on ne s’est pas borné aux chiffres qui pouvaient donner des 
chiffres exacts de la valeur delà racine cherchée et l’on a 

calculé celte quantité y avec toute l’approximation que l’on 
voulait obtenir pour la racine, et qui est ici celle à laquelle 
on peut atteindre par l’emploi des logarithmes. Si l’on s’é- 
tait arrêté à une valeur moins approchée de la racine de 
l’équation n-f- Az -+- llz* = o, il aurait fallu tenir compte 
dans le calcul de «'des termes n -+- A y H- 11'/* dont la somme 
n’eût pas été zéro. Mais, si l’on n’employait pas les loga- 
rithmes, ou si l’on voulait une approximation plus grande 
(pic celle qui aurait été d’abord obtenue parleur moyen, on 
devrait s’arrêter pour chaque correction aux chiffres dont 
l’exactitude pourrait être regardée comme certaine, afin de 
ne pas augmenter sans utilité la longueur des calculs , qui 
deviendraient alors fort pénibles. 


*) En changeant z' en — z ' el on remarquant que Ton ne devrait prendre 
que la plus petite racine de Toquai ion du second degré u ' — A' *'-f- IV ?'*== o, 
puisque l’autre racine serait plus grande que i, on aurait 


A' — y/A M — 4 IV u' _ 


>. IV 


A' -h v/A'*- 4»'u' 


, eu retranchant de cette quantité — » la différence 




..'(A' - y/À ,î -4BV) 4 BV* 

A'( \ B' u') A'( v/A’ s — 4 B'~m ’)’ 

H est facile de reconnaître , d’après les valeurs de A', B', u', que celle de la 
dernière quantité est moindre que 0,0000001 . 
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Méthode des approximations successives . 

21. Lorsque la correction z, qui doit être donnée par 

l'équation (6), est unequantité très-petite, si la valeur — 

qu’on obtient eu ne considérant que les deux premiers ter- 
mes, n’est pas assez approchée , on peut en avoir une plus 
approchée sans calculer une nouvelle équation. On a exac- 
tement, d’après l’équation (6), 



Soit « le quotient de la division de — u par A, ou une 
valeur approchée de ce quotient; en remplaçant z par a 
dans tous les termes du second membre, oir obtiendra une 
nouvelle valeur l> de z. On pourra alors remplacer z par b 
dans tous les termes du second membre, ce qui donnera 
une troisième valeur c de z ; ainsi de suite. 

22. Supposons, pour fixer les idées, que la quantités soit 
négative, et que tous les coeHicient-s A , B, C, D, E, etc., 
soient positifs. Dans ce cas, la valeur de zest moindre que 

— ^ • Donc, en remplaçant z par — ^ dans tous les termes 

du second membre, on aura une valeur trop petite de z. En 
remplaçant dans ces mêmes termes z par cette seconde 
valeur, on obtiendra une valeur trop grande; ainsi de 
suite. 

Four apprécier le. degré d’approximation des valeurs qui 
résulteront de ces corrections successives, désignons par P, 

O, K , etc., les coefficients -i-i-, etc., et soient acta -t-d 

A A A 

deux limites de la valeur dez. La valeur exacte sera moindre 
que celle qu’on obtiendra en remplaçant, dans le second 
membre de l’équation, z par a; et elle sera plus grande que 

a. 
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celle qu’on obtiendra en remplaçant z par « - 1 - d. Or, la 
difiérence des résultats de ces deux substitutions sera 

(2Pa+3Qa’-t-4jRa 3 + ...)£+(P+3Qa4-6Ra’-|- 

La quantité ô devant être trcs-pctite, et, en outre, d’un 
ordre de grandeur inférieur à celui de a, l’approximation 
de la valeur de z déduite de l’une ou de l’autre des deux li- 
mites a et a -(- ô dépendra principalement du premier 
terme de la difiérence des deux nouvelles limites; et elle 
sera sensiblement proportionnelle au produit ad. 

23. Si les coefficients des puissances de z dans le second 
membre de l’équation (ta) n’ont pas tous le même signe, on 
pourra obtenir deux limites de la valeur des, au moyen des 
deux limites a et a-f-d, en remplaçant z d’abord par a 
dans les termes positifs du second membre, et par a 4 - d 
dans les termes négatifs; et ensuite par a •+• d dans les 
termes positifs, etparadaœ les termes négatifs. Soient et ( z ) 
la somme des termes positifs, et — z ( 2 ) celle des termes 
négatifs; la différence des résultats des deux substitutions 
ou la limite de l’erreur des valeurs qu’on en conclura 
sera [ta 7 (a) z' (a) ] d -+- 7 [nr 77 (a) - 4 - n" (a)] d* -+- . . . . 

24. I.e mode d’approximation qui vient d’être indiqué 
peut aussi être employé dans des conditions différentes de 
celles quenous avons supposées , comme le montre l’exemple 
qui suit. 

L’équation par laquelle on détermine le diamètre d’une 
conduite d’eau cylindrique, est la suivante : 

LO 0 ’ LQ’ 

D s — 0,000. 095.94 — o >082.6^ D — 0,002. 22 — = 0. 

L est la longueur de la conduite supposée rectiligne et d’un 
diamètre uniforme; Q est la dépense, ou le volume d’eau 
qui s’écoule en une seconde; H est la pression à l’orifice de 
sortie , cette pression étant représentée par la hauteur de la 
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colonne d'eau propre à la produire; D est le diamètre in- 
connu. Les coefficients numériques ont été obtenus en sup- 
posant les longueurs rapportées au mètre pour unité. 

Cette équation a une seule racine positive qui est celle 
que l’on doit se proposer d'obtenir. Or, les coefficients du 
deuxième et du troisième terme étant très-petits, on peut 
négliger ces termes pour une première approximation. On 
a ainsi une valeur trop faible, qui est 

• D = ^0,002.22^- = 0,295 

Soient L = 757"', Q = o mmm ,o 8 g, H = 1"'. L’équation de- 
vient par ces valeurs , 

D 5 —0,006.464 D* — o, 000. 654 D — o,oi3.3i =0. 

Comme la question n exige pas une approximation de plus 
de quatre chiffres, on pourra prendre les logarithmes avec 
cinq chiffres décimaux seulement. 

En mettant pour D, dans le deuxième terme et dans le 

b - ' 

troisième, la valeur 0>,oi3.3i, on trouve : 


logD = i , 624.84 


logo, 000 . 654 = 4,8i5.58 

4 , 440.42 

% 

O 

O 

O 

log D 2 = 1 , 249.68 

„ . 

logo, 006. 464 = 3,8io. 5o 

3,o6o. 18 

o,ooi . i4g 
o,oi3.3i 

log de la somme 2 , i 68 . 35 

5 

0,014.735 

V 7,633.67 



L’excès du logarithme de la nouvelle valeur de D sur ce- 
lui de la valeur précédente étant 0,008. 83, 011 ajoute celte 
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quantité au logarithme précédent de ( 0,000. 654 ) I), et le 
double de cette quantité au logarithme de (0,006.464) D’- 
On a alors 

r , 

log (0,000.654) U — 4>449' 2 ^ 0,000.281 

log (0,006.464) D 2 =3,077.84 0,001.196 

0,0 i3. 3i 

log de la somme 2,169.88 0,014.767 

\J 1 ,633.98 

L’excès de la nouvelle valeur de log D sur la valeur pré- 
cédente est o,ooo. 3 1, d’où l’on conclut les nombres ci- 
après : 

4 , 44 g . 56 0,000.281 

3,078,46 0,001.198 

o,oi 3 . 3 i 

2,169.94 0,014.789 

1 ,633 99 

Le nombre correspondant au dernier logarithme 11e 
diilère pas, dans les quatre premiers ch i lire s , de celui qui 
était donné par le logarithme précédent; et il est o, 43 o 5 . 
Si l’on essaye, d’après cela, pour D, le nombre o, 43 o 6 , 
on trouvera que le logarithme de ce nombre surpasse 
1,633.98 de 0,000.09 , et l’on aura pour cette valeur : 

5 log D = 2,170.35 
log(o,oo 6 . . .) D = 4 , 44 g . 65 
log ( o , 00 . . . ) D 5 = 3,078.64 


O 11 conclut de cette vérification- que le nombre o,43o6 
est un peu trop grand, et comme le nombre précédent 
o,43o5 était trop petit, la valeur de D est connue, à moins 
de o '", 0001 . 
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Exemples tic la résolution des équations transcendantes . 

« 

25 . Le problème connu eu astronomie sous le nom de 
problème de Kepler consiste dans la résolution de l’équa- 
tion suivante : 

ii — c sin u — Ç. 

Les quantités u, e et £ sont des angles, le nombre e est 
aussi la valeur de l’excentricité de l’orbite. L’angle £ est 
donné en parties de la circonférence, et le nombre e est 
donné en parties du rayon ; pour le convertir en parties de 
la circonférence , il faut le diviser par sin i" . 

On reconnaît , par la considération de la dérivée, que la 
quantité u — e sin u est toujours croissante avec //. Pour 
les valeurs o et i8o° de on a u = Ç. Quand £ est plus 
petit que i8o°, u est aussi plus petit que i8o H , et l’on a 
u>£. Quand Ç est plus grand que i8o°, u est aussi plus 
grand que i8o°, et on a u Ç. ■ 

20 . Soit e= o ,5 (ce cjui suppose que l’astre soit une 
comète) et Ç — 38 ° 27' 1 8", 7. L’équation sera 

o ,5 


— sin u — 38 " 27' 18", 7. 


o ,5 


Le log de est 5 , o 1 3 . 3 g 5 1 , qui correspond à un angle 

compris entre 28° et 29 0 ; et comme on doit avoir sin u > 7 , 
puisque m^>£, ljangle u doit surpasser £ de plus de 1 4°; 
en conséquence, je prends pour un premier essai a = 54 °. 
La différence u — ( 38 ° 27'i 8",7) est alors t5°32 , 4i ff ,3 ou 
55961 ", 3 ; et l’on trouve que le terme e sin u est 83435 ", 8 - 
Le nombre 54 ° est donc trop petit. Si l’on prend 
u — 38 ” 27' t8",7 -+- 83435", 8 ou u = 6t° 87' 54", 5 , le 
terme <?siu u est 90747", 4 5 cette seconde valeur de « est 
donc encore trop petite. On pourrait 1 augmenter de l’excès 
de 90747", 4 sur 83434^, 8 ; on calculerait alors la valeur 
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du terme e sin u ; il en résulterait une nouvelle augmen- 
tation de la valeur de u; et ainsi de suite. Mais, pour 
avancer plus rapidement, on calculera d’abord la différence 
u — £ par une proportion , comme il suit. On a trouvé qu’en 
supposant d’abord u — 55 g 6 i", 3 , la différence des deux 

quantités u — £ et esin u est 83435", 8 — 55 g 6 i ", 3 
où 2 j 474 ") 5 , et que, lorsqu’on augmente u — £ de la 
dernière quantité , la différence entre u — £ et e sin u est 
90-47^,4 — 83435", 8 ou 7311", 6; de sorte que cette diffé- 
rence a diminué de '}.’j 474 "'r> — 73xi", 6 ou 20162", 9. On 
supposera , d’après ces résultats , que , pour annuler la 
différence entre u — Ç et e sin u , il faut augmenter la pre- 
mière valeur de « — £ du quatrième terme de la proportion 

20 162", 9: 27474 ” >5 :: 27474 * ',5 : x , 

d’où 

x = 37437 ", 5 = to° *3^57 ", 5 et « = 64 ” 23' 57", 5 . 

Parvenu à ce point, on se servira de la méthode de 
Newton; en désignant alors par u' ce qu’il faudra ajouter 
à la dernière valeur de u , on aura 

, c sin u — ( u — Ç) e sin u — 93398", 8 

1 — e cos u 1 — e cos « 

Le nombre t, dans cette formule, doit être divisé par 
sin 1" ; il représente un angle de 206260", et, d’après la 
valeur précédente de m, on trouve u' — — 499 w - Cette cor- 
rection réduit la valeur de u — £ à 92899", 8 , dont le loga- 
rithme est 4î9b’8.oi48. La valeur correspondante de « 
est 64 ° i 5 ' 38 ", 5 , et le logarithme de e sin u est 4 ,968 . 01 35 . 
L’équation n’est pas encore vérifiée avec une exactitude 
suffisante; mais, si l’où diminuait l’angle >< de 1", le loga- 
rithme de u — Ç diminuerait de 47, et celui de e sin u di- 
minuerait seulement de 10; donc la différence de ces loga- 
rithmes diminuerait de 37501 , comme cette différence est 1 3 , 
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il faudra diminuer l'angle u de j 4 , ou o", 35 . On aura ainsi , 
finalement, u = 64° 1 5 ' 38 ", i 5 . 

Voici le tableau sommaire de la marche du calcul : 

pour e==o,5, log -j-î-j. — 5 ,oi3.395i. 

u = 54° 

“ — Ç = i5° 32'4 i",3=5%6i", 3 <?sin «==83435", 8 = a3° io'35",8 
« = 6 i° 37 ' 54 ", 5 e sin u = go 747",4 

difler. 7311", 6 

83435 ", 8 

55961", 3 21052747415 = 8,877.3596 

diff. 27474 ,5 cowpl-log aoi6a ,9 = 5,695.4470 

73 11 >6 somme 4>573.3oG6 i 87417 ", 5 

diff. 20162 ,9 i 011 > 0 ° 23 ' 57 ", 5 

u = G 4 “ 23 ' 57", 5 

u — Ç= 93398", 8 c sin «= 93007", 6 391", 2 

ecosa= 44563 ” \ 

—^, = 206265" 1 J6.702" 

sin 1" J 

log 3gi,2 = 2 , 592 . 3g88 

compl . log 161702 = 4 791 .284G 
compl. logsin 1" = 5 , 3 i 4 . 4 t 5 i • • 

somme _ ' 2,698.1085 499"= 8* 19" 
u = 64 ° i 5 * 33 ", 5 

» — ç= 92899", 8 

log («— Ç) = 4,968.0148 log c sin a= 4,968 oi 35 J-î=o ,35 

b = 64 °i 5'38 ", i 5 

27 . Soit l’équation 

c z — e x z= ax, 


. qui se présente dans la théorie de la chaînette; et pre- 
nons a = 7 , 384 , L'équation est vérifiée par x = o ; mais 
elle admet une autre solution positive , qui est celle qu’il 
s’agira de trouver. 

Le nombre e étant moindre que 3 , la valeur de x 
devra être plus grande que 2. Pour ,r = 3 , on trouve 

c x — 2 o,o 8553 , c~ x = 0,04978; 

d'où 

e x — c ~ x - — n-v — — 2, r 1 625. 
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Si l’on faisait x — 4 , le terme «' surpasserait 54 , et, comme 
le terme ax augmenterait seulement de la valeur de a , 
ou 7 , 384 , la différence e x — e~ x — 'ax serait positive. La 
valeur de x est donc comprise entre 3 et 4 - On voit, de 
plus, que cette valeur doit être beaucoup plus près de 3 que 
de 4 j c’est pourquoi on essayera 3 ,i. On trouvera alors 


e 1 = 22, 19795, e~* — o,o 45 o 5 , 

et la valeur de e x — e~ x surpassera la précédente de 2, 1 1715. 
Comme la quantité axaugmentera en même tempsde 0,7384, 
la différence e x — e~ x — ax sera eucore uégalive et 
égale à — 0,7376; mais, si l’accroissement de c x — e~ x — ax 
était proportionnel à celui de x, la valeur cherchée de x 
surpasserait peu 3 , 1 5 . 

Soite r — e~ x — ax= f (x) ,d’où/ v (x) =e x -he~ x — a 
et J" (x} = e x — e~ x . On a exactement 

/ (* + h) =/(.r) 4 - hf (.r ) 4- - /* (x + 0 I, ) ; 


donc , en posant./ (x -I- h) — o , 




/ 




h’ f" (x 
2 


Oh) 


/'(*) 


et, en négligeant le terme dépendant de h % l’erreur, pour 
les valeurs de x plus grandes que 3 , sera moindre que h*. 

Pour x = 3 ,i, on trouve, à l’aide des résultats précé- 
dents, h — 0,04 , donc x = 3 ,t 4 - En mettant cette valeur 

de x dans —-CLÛ, on obtient h'= 0,0079, d’où x= 3 , 1 479. 

f[*)‘ 

Cette valeur est un peu trop grande; mais, si on la dimi- 
nuait de o,oooz, elle deviendrait trop petite. Si I on vou- 
lait une valeur plus approchée, il faudrait calculer une 
troisième correction , et l’on obtiendrait, celte fois , la valeur 
de X avec huit décimales. 
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Du calcul des solutions réelles de. deux équations 
simultanées. 

28. Lorsqu'on doit calculer les solutions réelles de deux 
équations simultanées à deux inconnues, on peut construire 
les courbes qui résultent de ces équations considérées sépa- 
rément-, les couples de coordonnées des points d'intersec- 
tion des deux courbes seront les solutions demandées. Mais, 
si l’on veut une approximation plus grande que celle qui 
pourra être fournie par les opérations graphiques , on fera 
usage de la méthode suivante, analogue à celle de Newton 
pour les équations à une seule inconnue. 

Soient J'lx,y) = o et F (x, y) = o les deux équations 
proposées, et soient a et b des nombres peu différents 
de ceux qui forment une dés solutions réelles. Ou posera 
x = a u, y = b ■ +- e; les deux équations deviendront 

(,31 ( ^ a > b )+S'J n ’ b )‘ a +S'A a ’ *)•<’+ u =o, 

( F (a, b) , b). ti («, b\,p -t- V =-o. 

f ' x («, b) désigne la dérivée du polynôme f [x , y), prise 
par rapport à x, en considérant y comme une constante, 
et dans laquelle on a remplacé x par a eiy par b ; de môme 
J'y {a, b) est la dérivée prise par rapport à y , en regar- 
dant x comme une constante, et dans laquelle on a rem- 
placé x par a et y par b. U et V sont des polynômes dont 
tous les termes contiennent des puissances ou des produits 
des inconnues u et e, d un degré supérieur au premier par 
rapport à ces inconnues. 

Les valeurs qu'il s’agira d’obtenir pour u et v devant être 
très- petites, on négligera les quantités U et V; ou n’aura 
plus alors à résoudre que deux équations du premier degré. 
Soient h et k les valeurs qu’on déduira pour u et v de ces 
deux équations du premier degré, les valeurs correspon- 
dantes de x et de x seront a -+- h et b -y- k. Nommons a, et />, 
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ces valeurs; si l’on a lieu de craindre qu’elles ne soient pas 
suffisamment approchées , on calculera par leur moyen une 
nouvelle correction , de la même manière que la précédente. 
On continuera ainsi autant qu’il sera nécessaire. 

Les premiers membres des équations (i3) étant les résul- 
tats de la substitution de a -+- u et de b -f- v , pour x etj\ 
dans les premiers membres des équations proposées , les 
valeurs que prendront U et V, lorsqu’on y remplacera u et v 
par les valeurs approchées h et À", seront celles que prendront 
les polynômes f{x,j) et F ( x , y) quand on fera x—a t et y=b, 
Donc, si , en prenant successivement, dans les équations ( 1 3) , 
à la place de a et de b, d’abord a, et /> ,, puis les valeurs 
qu’on aura déduites de celles-ci, et ainsi de suite, les va- 
leurs des inconnues it et v déduites des deux équations du 
premier degré deviennent de plus en plus petites , ce qui fera 
converger les quantités U et V vers zéro ; les nombres a et b, 
a , et b ,, etc., convergeront vers un des couples de valeurs 
réelles de x et de y qui vérifieront les deux équations pro- 
posées. 

29. Soient les deux équations 

y 7 — 4 *’r+* t=o, 

y 7 -4 - xy — îi 1 — Gy 4- 5t — 3i3 = o. 

Sans construire exactement les hyperboles déterminées par 
ces deux équations, on voit qu’elles se coupent en trois 
points très-voisins de ceux qui résultent des intersections 
des asymptotes. On peut obtenir exactement les coordonnées 
des points d’intersections des asymptotes ; on a , pour un de 
ces points, x — 5, 25, 10 , 5 ; pour un autre, x— — 15, 

y = — 3o,et, pour le troisième, x = -h 5 , y — — io. 
Chacun de ces trois couples est une solution approchée des 
deux équations proposées. Considérons le premier. En rem- 
plaçant x par 5,a5 -t- n et y par io,5 -4- f, on obtient les 
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deux équations 

I - 4 - 21 (' — 42 « .+ P ’ — 4 « ’ = O , 

2 4 - 20,25 e + 4°>5 a-f-p’ + up — 2 « 2 = o ; 

cl en négligeant les termes du second degré , par rapport 
à « et e, on trouve 

21,75 124,5 

u = — — 5 p = — . 

1701 1701 

Si l’on veut une approximation plus grande, on mettra 
d’abord les valeurs de u et de e en décimales, pour la com- 
modité du calcul; en s’arrêtant aux millièmes, on a 

u = — o,oi 3 , p= — 0,073. 

Il faut remplacer, dans les équations précédentes, u par 
— 0,01 3 -H v par — 0,073 4 - e,, ce qui conduit à 

0,017653 4 - 20,854 ‘'i — 4 1 5896 “1 4 - pJ — 4 ^ = 0, 

0,001 ig -+- 20,091 p, 4 - 4 o ,479 «1 4 - p? 4 - «1 p, — 2 u * = o. 

On conclut de celles-ci , en négligeant les termes du second 
degré, et en se servant des logarithmes, 

«, = 4- 0.000.195.62, p, = — o,ooo. 453 . 44 - 

En se bornant à six décimales , 011 a 

x = 5 . 237. 196, ^=10,426.547. 

Pour les valeurs de u, et de , déduites des dernières 
équations , en y négligeant les termes du second degré , 
ces termes ne donneraient que des quantités moindres 
que 0,000.001. Or, si l’on avait les valeurs de x et y qui 
vérifient exactement les équations proposées, une altéra- 
tion d’une unité du sixième ordre décimal , sur l’une ou 
l’autre de ces valeurs, rendrait les premiers membres diffé- 
rents de zéro de plus de 0,000.001. On est donc en droit 
de regarder les valeurs de .r et de y comme exactes jus- 
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qu’au sixième chiffre décimal. On peut vérifier toutes les 
opérations en calculant, à l’aide des logarithmes , les ré- 
sultats de la substitution de ces valeurs dans les premiers 
membres des équations; on trouve qu’ils ne di fièrent de 
zéro que de quantités qui sont inappréciables par ce mode 
de calcul. 


Du calcul ries logarithmes et des quantités trigonomctriqnes 
par le moyen des différences. 

. 30. Le calcul des valeurs d’une fonction au moyen des 
différences s’applique avec avantage aux fonctions autres 
que les fonctions entières. 

Considérons, par exemple, la fonction J. (x);'on trouve, 
en formant les différences successives , 


1Li = L(x+/i) — L ( .r ) 

(h A * A- \ 

— Ml ; 4- v - — — . . . ) i 

y.r si' jr j 

I. (.r) = L (x + 2 A ) — aL(x + A) + \,x 


2 h 


2L| 

/ h 

H 

x 

)- 

1 H 

\ x 

{/r 

2 A 3 



b- 

x z 

-f- . 

■Y 


i ;; L x = L(x + 3 A) — 3 L(jr -f- 2 h) -f- 3L {x — (— A ) — Lx 
= 1 . (, + H)_,3I.(' + ^) + 3I.(, + î 


"(T? — -)': 


etc. 


On voit que, pour /t = i , et en prenant x — i oooo , 
ou x ioooo, les différences successives décroîtront très- 
rapidement, et celle du troisième ordre, d'après là valeur 
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* » 
du module M , sera moindre que 0,000000 000000 869. Il 

suit de là, que si l’on ne voulait avoir les logaritlimes 

qu’avec huit décimales, on pourrait négliger longtemps les 

différences du troisième ordre, et ne tenir compte que de 

celles du premier et du deuxième ordre. 

La formule (2) (page 4 ) permet d’apprécier l’erreur 
qu’occasionne, sur une valeur d’un rang quelconque , la 
suppression des différences d’un ordre donné. Dans l’exem- 
ple que nous considérons, en faisant n — 4« , et en calcu- 
lant, pour cette valeur, celle du terme 

n( n —i)(n— 2) 

1 . 2.3 

on voit qu'il n influera pas sur la huitième décimale du 
logarithme de 1 o4o ; il en serait à plus forte raison de même 
des différences des ordres supérieurs. On pourra donc cal- 
culer les logarithmes des nombres depuis 10 000 jusqu’à 
10040, en partant de celui de 10000 qui est 4 ? et en cal- 
culant les valeurs correspondantes de A I..r et A ! L.r, qui 
sont : 

à L x = o , 000 .043.427.27, 

- o, 000.000.004. 34 - 

II est nécessaire de prendre ccs valeurs avec quelques 
chiffres au delà de ceux que l’on veut avoir exactement dans 
les derniers résultats, afin d’éviter l’accumulation des er- 
reurs. En négligeant dans A* Lx les chiffres qui suivent le 
onzième, l’erreur qui en résulte sur le logarithme de 10040 

est celle du terme — A’ u pour n — 4o, et elle est in- 

férieure à une unité du huitième ordre. 

Lorsqu’on sera parvenu au logarithme de 1 040, il faudra 
le calculer directement afin d’apprécier le degré d’exactitude 
des opérations. On calculera aussi les différences première 
et seconde correspondantes, et l’on formera ensuite, par 
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leur moyen, les logarithmes des nombres consécutifs dans 

un nouvel intervalle 5 ainsi de suite. 

31 . Considérons aussi les différences de la fonction sin x ; 
011 a : 

A sinx = sin(x-|-/() — sinx = 2 sin 3- A cos (x 4- A h ), 

A 2 sin x = A sin (x 4- A) — A sin x = — (2 sin 7 A} 3 sin (x 4- A),. 
A 3 sinx = A ! sin (x-f- A) — A’sinx = — (2 sin - A) 3 sin (■*-+- A A), 

A* sinx = A 3 sin (x4- A) — A 3 sinx = (2 sin ~ A)* sin{x 4- 2 A), 

etc. 

Il est facile de reconnaître la loi de ces différences, et l’on 
voit qu’elles décroissent très-rapidement lorsque l’arc h est 
très-petit. 

On obtiendrait des formules analogues pour la fonction 
cos x. 

32 . On a 

A 3 sin x = sin (x 2 A) — sin (x -+- A ) — [sin (x 4- A) — sin x], 

et, en faisant x = (m — i)h, on obtient, par la valeur ci- 
dessus de A 8 sin x, la formule de Thomas Simpson pour le 
calcul des Tables de sinus, 

sin (m 4- 1) A — sin mh = sin (mh) — sin (m — 1) A 
— (2 sin \h) 7 sin mh. 

Si l’on pouvait avoir les valeurs exactes des différences 
A 8 sinx, et celle de A sin x pour x = o, on en déduirait les 
valeurs exactes des sinus de tous les arcs croissant depuis 
zéro par des intervalles égaux à h. Les erreurs des valeurs 
des sinus calculées successivement ne seront donc occasion- 
nées que par les erreurs de ces différences. 

Soit e l’erreur du premier terme de la suite des diffé- 
rences secondes ou — (2 sin 7 A)* sin h. Cette erreur pro- 
duira une erreur égale sur chacun des termes de la suite des 


Digitized by Google 



33 


COMPLÉMENT o’aLÇÈBRE. 


différences premières à partir du second, et ces erreurs des 
différences premières produiront, sur les termes de la suite 
des sinus, à partir du troisième terme ou sin 2 A, des er- 
reurs croissantes e, ae, 3e. 4 e , etc. 

Si le second terme de la suite des différences secondes est 
encore fautif de la même quantité e, il en résultera sur les 
sinus consécutifs, à partir de sin 3 h, des erreurs e, 2e, 
3e, ainsi de suite. On conclut de là qu’en supposant toutes 
les différences secondes fautives de la même quantité e, et 
que le dernier sinus calculé soit sin mh , l’erreur de ce sinus 
produite par les erreurs des différences secondes sera 


; ( m — 1 -4- m — 2 -+- m — % -+- . 


m(m — 1) 

. -+- 1 ) ou — 5 e. 


Quant à l’erreur qui résultera de celle du premier terme 
de la suite des différences premières, elle sera seulement 
multipliée par m. 

En général, lorsqu'on calcule une suite de quantités 
«o, m, , « 2 , etc., au moyen de leurs différences, si l’une des 
différences de l’ordre n , A" u p est fautive d’une quantité e, 
tontes les différences de l’ordre n — 1, depuis A" -1 , 

sont fautives de la même quantité , et les erreurs qui en 
résultent sur les différences des ordres inférieurs et sur 
les quantités qu’il s’agit d’obtenir, à partir de A n-! u p+t , 
A" -3 u p+3 ..., «,, + „ , croissent, pour chaque suite, comme les 
nombres figurés des ordres successifs. Si toutes les diffé- 
rences de l’ordre n depuis A " u p sont fautives de la même 
quantité e, l’erreur sur un terme quelconque delà suite 
«o , «i , u s , etc., d’un rang supérieur àp + n, par exemple 
«p+n+m-i, est la quantité e multipliée par la somme des m 
premiers nombres figurés de l’ordre n — 1 , ou le m ,e ' n ' nom- 
bre figuré de l’ordre n. ( V. le Traité élèm. d’ Alg. , p. 358.) 

Pour le calcul des sinus, lorsqu’on fait h — io ff , on peut 
obtenir les différences secondes 


— (2sin{/i) ! sin(.r -4 -h), 


3 
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avec seize décimales exactes, ou à moins de.-^- { (*) ; la 

différence première pour li— o est sin h , et elle est connue 
avec treize décimales. 11 suflitde calculer les sinus et cosinus 
jusqu’à 3 o° ou 108000", de sorte qu’on a, pour le dernier 
terme de chaque suite , m = 10 800. On conclut de là qu’en 
supposant que toutes les erreurs atteignent leur limite su- 
périeure, et qu’elles soient toutes dans le même sens, les 
erreurs des différences secondes n’occasionneront, sur au- 
cun terme de la suite des sinus, une erreur supérieure à 


10800 X 10799 


, ou o, 000, 000. oo5. 83 1.46. L’erreur du 


2 X 10” 

premier terme de la suite des différences premières, ou 
sin 10", produira sur le dernier sinus une erreur moindre 
que 0,000.000.001. 08. Il y aura donc au moins huit dé- 
cimales exactes dans toute la suite des résultats. 


(*) On a pour h = 10", 

sin j h < 0 , 000 . 024 - 2 '|0. 684. o 5 fv. . 


On conclu! de là 


i (î)' (0,o *°' o5), ‘ 

1 /h \ 3 

n 1 *J j < '° 9 000 * 000 • 000 • 0(H> °°3, 
si 11 J- h > 0,000 024 . 240,684 . © 5 a. 


1 . <0,000.048.481 . 368 . 1 10 
2 sin — h _ .• 

’â >0,000.048.481.368.104 

Ainsi la valeur de a ain J h est connue avec quatorze décimales exactes. En 
multipliant celte valeur par eïle-fiiôroc,, on aura dix-neuf décimales dans 
tous les produits partiels, cl Ton connaîtra le produit total avec dix-sept 
décimales. Le facteur sin r 4- h), par lequel le nombre constant (2sin£A)* 
doit être multiplié, étant toujours moindre que 1, la multiplication don- 
nera dix-huit chiffres à chaque produit partiel, cl le produit total n’aura 
pas moins de seize chiffres exacts Comme le nombre (2 sin | A) 4 a , d’ail- 
leurs, huit zéros entre la virgiHè et le premier chiffre significatif, l’erreur 
du facteur sin ( x -+• h) no pourra influer sur les seize premières décimales 
du produit, tant qu’elle n'aura pas dépassé elle-même une unité du hui- 
tième ordre décimal. 
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De l'interpolation. 

% 

33. On nomme interpolation , l’opération par laquelle, 
étant donné un certain nombre de valeurs d’une fonction, 
avec les valeurs de la variable auxquelles elles correspon- 
dent, on détermine les valeurs de cette fonction correspon- 
dantes à d’autres valeurs intermédiaires et données de la 
variable. 

Dans le cas où la fonction proposée est entière et de de- 
gré n , elle est complètement déterminée quand on connaît 
les valeurs qu’elle prend pour «-+- i valeurs différentes 
et données de la variable (n°7); et, lorsque ces valeurs 
de la variable sont en progression par différence, l’expres-- 
sion générale de la fonction est celle qui résulte de la for- . 
mule (3) [page 8]. Mais, le plus ordinairement, on a 
recours à l’interpolation pour calculer les valeurs d’une 
quantité dont on n’a pas 1 expression générale, et dont on 
connaît* seulement un certain nombre de valeurs particu- • 
lières. Le problème est alors tout à fait indéterminé, car 
il revient à trouver une courbe qui passe par des points 
désignés; et, quel que soit le nombre de ces points, il est 
toujours possible de les joindre par une infinité de lignes 
courbes qui diffèrent- les unes des autres d’une manière 
sensible dans l’intervalle de chaque point au suivant. 

Cependant, comme les fonctions dont les valeurs se suc- 
cèdent suivant une certaine loi, sont ordinairement sus- 
ceptibles d cire exprimées par des séries ordonnées suivant 
les' puissances de la variable, et qui sont convergentes dans 
un certain intervalle des valeurs de celle variable, on con- 
çoit qu’on peut les assimiler aux fonctions entières, et 
obtenir ainsi des résultats suffisamment approchés. De cette 
manière, on choisit, pour la courbe dont l’ordonnée re- 
présente la fonction que l’on a à considérer, une courbe 
parabolique, et qui est la plus simple parmi toutes celles 

3. 

'J/ , ' ■ 

‘ t ' 4 * « , 
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36 

qui peuvent unir les points donnés. Mais , pour que oc 
mode d’interpolation conduise à des valeurs qui ne soient 
pas trop inexactes, il faut n’en faire usage que dans des 
intervalles assez, resserrés; et, quand on emploie la for- 
mule (3), il est nécessaire que les différences successives 
décroissent assez rapidement pour que l’on puisse négliger 
tomes celles qui dépassent un certain ordre. 

34. En faisant, dans la formule (3), X — x 0 = h', 
u x — «o = A' « , et en écrivant simplement u au lieu de u 0 , . 

il vient 


(>4) 



h' h' œ — i 

— bu -l — 

h 1.2,/j- 


A’ u 


— u -h . . . 

I.3.3.È 1 


Cette formule fera connaître la différence A'm, pour un 
intervalle h ', au moyen des différences relatives à l’inter- 
valle h. 

Si les différences seconde et suivantes sont assez petites 
pour qu’on puisse les négliger, la différence A' u sera pro- 
portionnelle à l’intervalle h',. C’est ce qui a lieu , par 
exemple, pour les usages les plus ordinaires des Tables de 
logarithmes. Mais il se rencontre aussi des cas, relative- 
ment à ces Tables, où l’on ne peut pas négliger tous les 
termes de la valeur de A'u qui suivent le premier, comme 
on le voit par l’application suivante, qui est indiquée par 
l.acroix. 

35. Supposons qu’on veuille obtenir le logarithme de 
3, 1 4 1 - 5t)2.653 .6, par le moyen d’une Table contenant les 
logarithmes des nombres depuis t jusqu’à iooo, avec dix 
décimales. On regardera les logarithmes contenus dans 
cette Table , comme les valeurs données de la fonction u x , 
les nombres comme celles de X , et l’on formera le tableau 


suivant : 
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SOUBRE3 . 

LOGARITHMES. 

DIFFÉRENCES l r ". 

IMITÉ R. 2'*. 

DIFF. S*‘. 

DU, 

ï l/â 

H 

o. 49 6 '92<J.648 .i 

H-I .380.905.7 

-4.376.9 

-1-27.7 

.* _ . 

-3 

3, i5 

.o, 498.3 io. 553.8 

H-l .376.528.8 

—4.349.2 

+27.4 


3,iG 

0,499.687.080.6 

-l-I .372.179.6 

-4.321.8 



3,-7 

o,5oi .059. 262. a 

-H 1 . 367 . 857 . 8 

.. 



j 3 ,i 8 

OjSoo. 427. 120.0 






En prenant quelques logarithmes consécutifs de plus, on 
trouverait encore — 3 pour la différence quatrième; de 
sorte que, lorsqu’on s’arrête à dix chiffres décimaux, les 
différences cinquième et suivantes sont nulles , et l’expres- 
sion de A' u, pour le nombre proposé, doit se terminer au 
quatrième terme. 

On a , d’après le tableau ci-dessus : 
u = 0,496.929.648.1, 

lu = ■+■ 0,001 . 380.905.7, y u = — o ,000.004 . 376 g, 
A 1 u = -f- o , 000 . 000 .027.7, - y u = — o , 000 . 000 . 000 . 3 ; 


et comme 

/i = o,oi, h ' = 0,001 .592.653.6, 
on obtiendra 



o , 1 59 . 265 . 36 , 


h'—h_U' 1 

2/1 2 h 2 


0,420.367.32, 


h ' — 2 h h' 2 

3 h = 3 Â“ 3 ~ 
h ' — 3 /i h ' 3 

4// 


— 0,613.578.21 , 


— 0,710. i 83 . 66 . 


Avec ces valeurs, il sera très-facile de mettre en nombre la 
formule (i 4 ), qui donnera A' u = 0,000. 220.224 • 5 , et, 
par conséquent , 

log 3 , i 4 t . 5 g 2 . 653 .6 = 0,497 .149.872.6. 
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36. Ce que l’on vient de dire de la formule (i4), par rap- 
port aux logarithmes des nombres, s’applique aussi aux lo- 
garithmes trigonomélriques , avec cette seule différence, 
que les quantités h et h' sont alors les dilférences des arcs , et 
non pas celles des nombres auxquels correspondent les lo- 
garithmes, et qui ne sont pas connus. Mais on peut remar- 
quer que les quantités h et h' n’entrent dans la formule (i 4 ) 
que par leur rapport, et que, pour des dilférences très-pe- 
tites, le rapport des dilférences des arcs ne diffère pas sen- 
siblement de celui des différences des quantités trigonomé- 
triques. On a, en effet, en considérant, par exemple, les 
sinus, 

sin(rt-4-z) — sin a sin-~zros(/ï -J-'-j*) 
sin;u-+-z') — sin« sin ÿz' cos (a -+- jz'J 

Si les arcs z et z' sont très-petits, et si, de plus, l’arc a 

i , . . , „ , cos (a -t- 4- z) . 

n est pas trcs-voisin de qo ”, le rapport } — — est tres- 

. r y n cos(<7 4-iz') 

peu différent de l’unité (*) ; et , en même temps, le rapport 
S1 " ; Z . est sensiblement égal à 

SID y Z " z' 

Lorsque les arcs a -f- j z et a -f- \z' sont très-voisins de 
90 °, le rapport de leurs cosinus peut être très-différent de 
l’unité, et il n’est plus permis de substituer le rapport des 
dilférences des arcs à celui des différences des sinus; mais, 
dans ce cas , les petites différences des arcs ne font varier les 
• sinus et leurs logarithmes que dé, quantités insensibles au 
degré d’approximation des Tables ; ainsi il n’y a pas lieu de 
calculer ces quantités. 

37. Quand les arcs sont très-petits , on 11 ’aurait pas une 
exactitude suffisante pour les logarithmes deg sinus ou des 
tangentes, en se bornant au premier terme -de la valeur de 

(*J. En cherchant dans les Tables la différence des logarithmes des cosinus 
tics arcs 89° cl 8j>° o' 10", on voit que le rapport do ces cosinus est moindre 
que t ,oo 3 . 
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A'», parce que les différences premières des logarithmes 
varient alors très-rapidement. Mais, au lieu de se servir de 
la formule (i4) , il est préférable d’employer le moyen qui 
est indiqué dans l’avertissement placé en tète des Tables de 
Callet, et que je ferai connaître ici, par un exemple ; pour 
la commodité du lecteur. . 

Supposons qu’on veuille avojr le logarithme du sinus de 
a3.", 63. Les arcs très-petits étant sensiblement proportion- 
nels à leurs sinus, on aura , à très-peu près, 

sin a3",63 = sin 23" X 

20 

On prendra donc, dans les Tables qui donnent les loga- 
rithmes des sinus et des tangentes des arcs de seconde en 
seconde, le logarithme du sinus de 23"; on y ajoutera lç 
logarithme de 23,63 et le complément du logarithme de 23 ; 
la somme diminuée de to sera le logarithme demandé. On 
suivra uite marche analogue pour trouver l’arc au moyen 
du log sin , quand cet arc sera très-petit. 

38. La formule ( 1 4 ) s’applique seulement au cas de l’in- 
terpolation d’après des valeurs connues , correspondantes à 
des valeurs équidistantes de la variable. Dans le cas plus 
général où l’on connaît n-h i valeurs d’une fonction, cor- 
respondantes à un égal nombre de valeurs quelconques don- 
nées de la variable , en supposant toujours la fonction en- 
tière et du degré n, on peut poser 

Uj — A + Bar + Car 3 -4- . . . -4- K.x"; 

on a alors , pour déterminer les coefficients A , B, C , . . . , K , 
les n -4- i équations suivantes : 

A -4- Bx ft -f- Car \ -4- ...-(- K-X^', 

H i — A -4- Bx, -|- Cx| -4- . . . “4" Kx" , 

■ K, — A -4- Bx, -4- Cxj - 4 - . . . -4- K x" , 

u n zz: A -4- IJéffl -f- -f- . . . -f- K. Jc n * 
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D’après la loi des valeurs des inconnues dans les équa- 
tions du premier degré, celles qu’on obtiendrait pour les 
coefficients A, B, C,..., R, ne contiendraient les quantités 
«o, u n qu’au premier degré et dans les numéra- 

teurs seulement, de sorte qu’il en résulterait une expression 
de cette forme , 

U-X Xg U u —J- X| li y -f- X, Uj . . . -f- X„ Mfl \ 

les fonctions X 0 , X,, X s ,..., X„ ne dépendant que de'la 
variable x et de ses valeurs x 0 , x, , x, x 0 . 

Cette expression de u x donnera n r = u„ pour x = x 0 , si 
l’on a, par cette valeur de x, 

X.( — i , X, — o, X, — o, . . . , X, — o j 

de même x = x, donnera u x — u t , si l’on a alors 

X, — o , X, — i , X, — o , . . . , Xjt — o y 

x = X| donnera u x — u,, s’il en résulte 



enfin x — x„ donnera «,= «„ , si l’on a par cette valeur 
de x, 

X — o, X i “ o, X, ~ o, . . . , - X„ ~ i . 

En considérant par colonnes le tableau de ces conditions, 
on voit d’abord que la fonction X 0 doit devenir nulle lors- 
qu’on donne à x toutes les valeurs comprises dans la série 
x 0 , x, , x,,..., x„, excepté la première. Or on satisfait à 
cette condition de la manière la plus simple , en prenant le 
produit 

(x — Jf, )(* J,) . . . (x — x„). 

Lorsqu’on fait x = x 0 , ce produit devient 

( Xÿ i J (x, Xi J . . . Xn ) j 

donc, en divisant le premier produit par le second, le 
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quotient, qui deviendra l'unité pour x — x„ , remplira 
toutes les conditions relatives à la fonction X 0 . On pourra 
donc prendre 

x = {■*■ — *■)(* — *0 • • • (* — x„) 

" (x, — x,)(x, — x,) . . .-(X, — x„) 

On trouvera de même 

x _ (x — x,)(x — X») ... . (x — X„) 

‘ — (x, — x, )(•*,-— x,) .. . (x, — x„) 

X — ( J — '*»)(* — X,) . . . (x — x„) 

’ 2— (Xj — x,)(x, X|) . . . (x, — x„) 

x _ (x — x„)(x — X,) ■ ■ ■ (x — x„_, ) 

(x„ X„){ X„ ) > « • ( X n X„_i ) 

On aura ainsi la formule suivante : 

(x — x,)(x — x,)...(x — x„) (x — x„)(x — x,)...(x — x„) 

ll x — rr r — ; : u c “T" 7 77 ; — 77 \ u i 

(■X 0 — X, XjJ. « ,(x^~-X n ) . .(«ï| X n ) 

(x — X.) (x -r- X,)...(x — x„_,) 

. (X„ — X 0 )(X„ — X,)...(x„ — X„_,) 

r 

Cette formule, duo à Lagrange, et remarquable par son 
élégance, est surtout très-commode pour les applications, 
puisqu’on peut calculer immédiatement tous les termes par 
les logarithmes. 

39. 11 est à remarquer que les considérations par lesquelles 
on a formé les valeurs des quantités X 0 , X, , X s , . . , , X„ 
sont bien propres à montrer comment l’interpolation est un 
problème indéterminé , lorsqu’on ignore la forme de la 
fonction d'où dérive la série des nombres donnés. En ellct, 
le produit des binômes x — x — x t ,..., x — x n n’est pas 

la seule fonction susceptible de s’évanouir pour les valeurs 
x = x t ,= X t ,...,= x n . Si l’on sort des fonctions algé- 
briques , on trouve d’abord l’expression très-simple 
sin /> (x — x, ) sin y (x — x, ) . . . sln / (x — x„ ), 


Digitized by Google 



r.OMPLKMKIVT d'ALGKBKF. 


qui jouit de la même propriété, quels que soient les nom- 
bres p, î; on pourra donc poser encore 



sin p (x — x, ) sin 17 (x — x,) . . . sin t(x — x„ ) 
sin p (x 0 — x.) sin q (x 0 — x,).. . . sinr(x„ — x„) ’ 


sin p' {x — x„) sin q’ (x — x,) . . . sin f'{x — x„) 
sin p' (x, — x„) sin </' (x, — x,). . .sin f' (x, — x„) ’ 


etc., 


en observant que les nombres/»', <7',..., t' peuvent être dif- 
férents des nombres/», q,..., et ainsi de suite pour toutes 
les autres fonctions X. 

Si de pareilles expressions s’accordent avec celles qui 
précèdent pour les valeurs de x comprises dans la série x 0 , 
x ,, x,,..., X n , elles en diffèrent beaucoup pour les valeurs 
intermédiaires dès que les arcs ne sont plus assez petits pour 
être sensiblement proportionnels à leurs sinus. 
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NOTE I. 


Sur la résolution do l'équation ax' l -\- bx -t- e = o , lorsque 
le produit ac est un nombre très-petit. 

•Lorsqu’on a à calculer les racines d'une équation du 
second degré, si le produit ac est un nombre très-petit, la 
quantité v o 5 — 4 nc est très-peu différente de />,■ cl, pour 
avoir avec une approximation suffisante la racine 

• , — b y b 1 — Aac 

x = — , 

2 a 

il est nécessaire de pousser l’extraction de la racine carrée 
de la quantité b * — 4 ac jusqu’à une approximation plus 
grande que celle qui peut être obtenue par l’emploi des 
logarithmes; mais, dans ce cas, ou peut faire usage de la 
série suivante : 

c ac 3 2 a 1 c 3 

b b' b 3 

On obtient cette série par la formule du binôme, pour le 
cas de l’exposant fractionnaire, et en développant, suivant 

cette formule, la quantité (&’ — 4 ac)*. Mais, le plus or- 
dinairement, il suffit de prendre le* premier terme ou les 
deux premiers termes de la série , et on les obtient aisément 
par les procédés élémentaires. 

En multipliant les deux termes de l’expression de x 1 par 
\Jb * — 4 ac ■+• b, 011 a 



V b 1 — 4 UC +■ b 


la quantité \jb l — 4 «c étant peu différente de b , si on U 
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remplace par b, ou aura une valeur approchée Je x', et celle 

•valeur sera — T • 
b 


L’erreur de cette valeur approchée sera 


•?c 


Or, cette expression est équivalente à 

c {\Jb * — 4 ac — b) — 4 ac " 


v^’-4 


ac- 


r * 


b^b* — 4 ac + b) b{\'b‘ — 4 ac - 4 - b) 

eu remplaçant sjlr — 4 ac par b dans le dénominateur de la 
dernière expression, on aura, pour valeur approchée de 


ne’ , , c ac' 

F’ donc x =~b~-F' 


cette expression, 

On pourrait trouver, d’une manière aualogue, les termes 
suivants de la série 5 mais les transformations seraient plus 
compliquées. 


La racine 


b — ^ b 1 — ■ 4 ac 


pourra se déduire de la ra- 


cine x', puisqu’on devra avoir x' x" = — L’erreur re- 
lative de la racine x " , calculée d’après cette relation, sera 
la même que celle de la valeur approchée de x'. 
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NOTE II. 


Surfes approximations numériques , et. sur une abréviation 
fie fa division. 

Lorsqu’on veut apprécier l’approximation d’un résultat 
numérique qui n’a pu être calculé avec une entière exac- 
titude, on doit moins considérer la valeur absolue de l’er- 
reur, que sa valeur relative; c’est-à-dire le rapport de 
l’erreur au nombre qu’on avait à calculer- 

Désignons par A un nombre inconnu, et soit À' une 
valeur approchée de ce nombre; l’erreur absolue sera la 

A A' 

dilTérence At— A', et l’erreur relative sera le rapport — • 

Nommons a ce rapport, qui devra être une très-petite frac- 
tion de l’unité ; on aura A'= A (t — a). Si le nombre A' 
a été calculé avec m chiffres à partir du chiffre des plus 
hautes unités, de telle sorte que l’erreur absolue A — À' soit 

4 X A' 

moindre qu’une uni té du m i, T r chiffre, l’erreur relative — — — 

sera moiftdrc que — ~ ; mais, pour que l’on soit assuré que 
l’erreur absolue n’est pas d’une unité du m Umt chiffre , il ne 
’ suffira pas que l’erreur relative soit moindre que et 

elle devra être moindre que (*). 


(*) l/erreur relative de la valeur approchée A' no changera pas si celle 
valeur et le nombre A sont rendus tous deux plus grands ou plus petits 
dans un môme rapport. On pourra donc supposer que le nombre A' de m 
chiffres est un nombre entier, la partie entière de A ayant aussi m chiffres. 

. j .A — A' t • * ‘ . 

Alors, si — t < ■ 1 on pourra avoir A — A >i: mais si 

A i o m— 1 1 

A — A' t 

— ; — < 1 on aura nécessairement A — A' < i. 

A io*" 
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Si A' et B' sont les valeurs approchées de deux nombres 
A et B dont on doit calculer le produit, et si les erreurs 
relatives de ces valeurs approchées sont a et S, celle du 
produit A' B' sera 1 — (1 — «) (1 — 0), ou a -t- 6 — a 6 ; si 
les valeurs A' et B' étaient toutes deux approchées par ex- 
cès , l’erreur relative du produit A' B' serait a 4- S -4- a ë ; 
mais la quantité a o étant très-petite par rapport à a et à S , 
qui doivent être de petites fractions de l’unité , on pourra 
admettre, dans tous les cas, que la limite de l’erreur rela- 
tive du produit est la somme des limites des erreurs des 
facteurs; et, quand celles-ci ne seront pas du même ordre, 
on poivra 11e tenir compte que de la plus grande (*). 


(*) Supposons les quantités a et 6 moindres que 0,001 ; on sera assure 
que l'erreur relative du produit sera moindre que 0,002001. Serait-on 
fondé à prétendre qu'elle pourra être plus grande que o,oo 2 7 surtout si les f 
quantités a et 6, comme cela arrivera toujours , sont non pas t-gales à o,oüi, 
mais toutes deux moindres que cette fraction. Il doit paraître .bien clair 
que, si l'on voulait pousser lu rigueur du raisonnement jusqu'à tenir 
compte du terme oc€ dans l'estimation de l'erreur du produit, on ne ferait 
que se créer, en pure perte, des embarras et des diüicullés pour la pra- 
tique. Mais, afin de montrer d'une manière frappante combien il serait peu 
judicieux de vouloir s'astreindre à cct e gard à une rigueur absolue, suppo- 
sons que les parties entières des nombres A et U aient quatre chiffres, et 
qu'on ait calculé ces nombres à moins d'une unité, pour que çe <0,001 et 
£< 0,001. Supposons, de plus, les valeurs A' et B' approchées par excès. 

Les erreurs a et 6 seront respectivement moindres que et donc , pour 

que l'erreur relative a-Hê-f-aS du produit fiU plus grande que 0,002. il. 
faudrait que l'on eût à fortiori 


Î + v,0 " a ' 


on devrait avoir aussi 


H > 10’, d'où < -- 
1} io‘ 


-> o,ooi', • d'où A < 1 1 xi 1 . 


A 1000 A 

J‘ar une raison semblable, on devrait avoir egalement Ii <1001. Ces coudi- 
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Si l'on doit calculer le quotient des deux nombres A cl R, 
dont les valeurs approchées sont A' et R', l’erreur relative 

du quotient -, ou sera x — ou — g- SoH y 

une fraction plus grande que chacune des deux fractions « 
et £. Si les deux nombres A' et R' sont approchés de A et 


de R dans le meme sens , la fraction 


a — 6 


sera moindre 


que — g ; et si les nombres A' et R' sont approchés de A 


tions seules ne suffiraient pas, et il faudrait, en outre , que les différences 
A' — A et B' — B fussent toutes deux très-voisines de l imité. 

Généralement, si A' a m chiffres et B', n chiffres, pour que l’erreur rela- 
tive du produit A' IV fût plus grande que — ^ — - H — ^ f i il faudrait que 

Ton eût 


> 


A B AB io'« 


et de l?i , à cause de B > 10 " -• , d'où — < , on conclurait 

B io n -‘ 


4 (ih — |> — - — -i A < io m -* 

A \ io"-' / io*-' 


io™ 
*o n J 


on devrait avoir pareillement 


D ^ n 10 

R < ÏO" H 

ro m 


Les illégalités 

x -t- 6 ■+■ w6 


-s a < 


l 0 m— i lo' 1- ‘ I0 m_l 

exigeraient de plus les conditions suivantes : 


ê< 


(, , 1 1 

\ "> * 1 « 6 1 


, *> 1 • 

\ 1 0" ‘ j 

’ io m —' ' 

{ m— V 

^ IO"-' *■ 


et puisque 

A r — A B' — B ' 

a = — r — » 6 = — ^ — » A > i o'"~ • , B >io"-»*^ 


on conclurait 


■ A > 


io w -' 


10"“' -t- i 

B # - B > 


ou A' — A > i — 


,o«-* i 
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4 » 

et de B en sens contraire, la fraction \ sera moindre 

9 I — 6 

O fy 

que ■ ■ ■ ■ D’ailleurs, S élant une fraction très-petite, on 

peut substituer aux limites — — s et > les limites y 

et iy {*). 

Supposons qu’on veuille obtenir exactement les m pre- 
* . A 

miers chiffres du quotient — ) à partir du chiffre des plus 

hautes 'unités, l’erreur relative du quotient approché devra 

être moindre que P ar conséquent, il faudra que les 

erreurs a et 6 soient aussi moindres que cette fraction si 
elles sont dans le même sens, et moindres que la moitié 
de la même fraction si elles sont de sens contraire. La 
première condition sera remplie, si les nombres A' et B' 
sont calculés avec m + 1 chiffres (**)•, et, pour la seconde, 

(*) Soit y = ce qui suppose que les nombres A' et B' aient été 

calculés avec m~\- i chiffres, el de manière que Terreur ne soit pas d'une 
unité de m •+• l^ème chiffre ; pour que Ton eût — H faudra il que 

i i 

oc > — — — ; 

io m IO*-’" 

et de là on conclurnil, en supposant, comme précédemment , que le 
m-f- i*«a* chiffre de A exprime des unités entières, 

i i * 

A<io* + H ■> A — A' > i 

IO w — I IO* 

(••) Si le premier chiffre du quotient doit ôtre moindre à la fois que 
le premier chiffre de A et que le premier chiffre de B, il suffira de calcu- 
ler A' et B' avec m chiffres. Et» effet, soient a et b les premiers chiffres de A 
et de B; les erreurs a et 6, en calculant A' et B' avec m chiffres , seront res- 
pectivement moindres que el - — ;• Donc , en supposant b t 

1 it.i «»"■“* 
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il suffira , dans tous les cas, que ces nombres soient calculés 
• avec tri -H a chiffres (*)■ 

Occupons-nous, maintenant, de la division des nombres 
A' et B', calculés comme nous venons de le dire-, si l’on 
donnait ces nombres, ou seulement l’un d’eux, avec plus 
de chiffres qu’il n’est prescrit dans ce qui précède , on 
pourrait supprimer les chiffres excédants. 

Quand on a obtenu le premier chiffre du quotient, puis- 
qu’il n’en reste plus à calculer que m — i , on peut con- 
server un chiffre de moins au diviseur, en 11e prenant aucun 
chiffre nouveau au dividende; pareillement, quand on a 
. trouvé le second chiffre du quotient, on peut supprimer 
uu second chiffre à la droite du diviseur; ainsi de suite. A 
la vérité , on commet ainsi des erreurs successives qui 
s’ajoutent, et comme chacune d’elles a seulement pour li- 
mite l’unité de chiffre du quotient, l’erreur totale 

peut être de plusieurs unités de cet ordre; mais, en ne 
commençant à supprimer un chiffre à la droite du diviseur 
qu’ après avoir obtenu les deux premiers chiffres du quo- 
tient, chacune des erreurs produites par les altérations 
successives du diviseur sera moindre que l’unité du 
m -+- i' emr chiffre du quotient; et, si l’on n’a pas à calculer 
en tout plus de douze chiffres, l’erreur finale sera moindre 
qu’une unité du dernier chiffre. Cette erreur augmentera 
le quotient, puisque le diviseur sera successivement dimi- 
nué; l’erreur qui résultera de ce qu’on négligera tous les 



i 


* 

i 


l’erreur absolue du quotient est moindre que «t puisque le pre- 

mier chiffre de t] est moindre que b , la fraction ^ est moindre qu’une 
unité de m iima chiffre de 7. 

(*) Toutes les fois que le premier chiffre de A est autre que * , en cal- 
culant A' avec m -+- 1. chiffres , on a v < : — 

• 2 . 10 "* 


4 
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chifl’res du quotient au delà du m“ mr , diminuera, au con- 
traire, le quotient; ainsi, ces deux erreurs se compense- * 
ront en partie (*). 

A' 

(*) S’il fallait tenir compte de ce que le quotient — pourra différer du 
À 

quotient de près d’une unité du m iime chiffre; comme on ne serait pas 

D • 

assure du sens de Perreur à Pétard du quotient des valeurs approchées , on 

ne pourrait , pas répondre que Perreur, par rapport au quotient — , ne 

s'élèvera pas à deux unités du m iim * chiffre. Pour éviter cette incertitude, il 
faudrait calculer les valeurs approchées et leur quotient avec un chiffre de « 

plus. Mais H sera ordinairement d'une très-médiocre utilité de s'astreindre 
à une exactitude aussi minutieuse. . 


FIN. 
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